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Abstract 

Radiatively-driven transfer flow perpendicular to a luminous disk is examined in the relativistic regime 
of {v/cY, taking into account the gravity of the central object. The flow is assumed to be vertical, and 
the gas pressure as well as the magnetic field are ignored. Using a velocity-dependent variable Eddington 
factor, we can solve the rigorous equations of the relativistic radiative flow accelerated up to the relativistic 
speed. For sufficiently luminous cases, the flow resembles the case without gravity. For less-luminous or 
small initial radius cases, however, the flow velocity decreases due to gravity. Application to a supercritical 
accretion disk with mass loss is briefly discussed. 
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1. Introduction 

Mass outflow from a luminous disk is a clue to the for- 
mation mechanism of astrophysical jets and winds in the 
active objects. In particular, in a supercritical accretion 
disk, the disk local luminosity exceeds the Eddington one, 
and mass loss from a disk surface driven by radiation pres- 
sure would take place (see Kato et al. 1998 for a review 
of accretion disks). 

So far, radiatively driven outflows from a luminous 
disk have been extensively studied by many researchers 
(Bisnovatyi-Kogan, Blinnikov 1977; Katz 1980; Icke 1980 
Meha, Konigl 1989; Misra, Meha 1993; Tajima, Fukue 
1996, 1998; Watarai, Fukue 1999; Hirai, Fukue 2001 
Fukue et al. 2001; Orihara, Fukue 2003), and by numer- 
ical simulations (Eggum et al. 1985, 1988; Okuda 2002 
Ohsuga et al. 2005; Ohsuga 2006). In almost all of these 
studies, however, the luminous disk was treated as an ex- 
ternal radiation source, and the radiation transfer in the 
flow was not solved. 

Radiation transfer in the disk, on the other hand, was 
investigated in relation to the structure of a static disk 
atmosphere and the spectral energy distribution from 
the disk surface (e.g., Meyer, Meyer-Hofmeister 1982; 
Cannizzo, Wheeler 1984; Shaviv, Wehrse 1986; Adam et 
al. 1988; Hubeny 1990; Ross et al. 1992; Artemova et al. 
1996; Hubeny, Hubeny 1997, 1998; Hubeny et al. 2000, 
2001; Davis et al. 2005; Hui et al. 2005). In these studies, 
however, the vertical movement and the mass loss were 
not considered. 

Recently, mass outflow as well as radiation transfer have 
been examined for the first time in the subrelativistic 
(Fukue 2005a, 2006a) and fully relativistic cases (Fukue 
2005b). In the latter case, it is pointed out some sin- 
gular behavior inherent in the relativistic radiative flow 
(e.g., TuroUa, Nobih 1988; TuroUa et al. 1995; DuUemond 
1999). When the gaseous flow is radiatively accelerated 
up to the relativistic regime, the velocity gradient becomes 



very large in the direction of the flow. As a result, the ra- 
diative diffusion may become anisotropic in the comoving 
frame of the gas. Hence, in a flow that is accelerated from 
subrelativistic to relativistic regimes, the Eddington fac- 
tor should be different from 1/3 even in the optically thick 
diffusion limit. 

In order to avoid the singular behavior of such a rel- 
ativistic regime, for a plane-parallel case, Fukue (2006b) 
proposed a velocity-dependent Eddington factor, which de- 
pends on the flow velocity v: 

1 + 2/3 



/(/?) 



(1) 



where (3 = v/c. In Fukue (2006b) this form (1) was 
adopted as the simplest one among various forms, which 
satisfy several necessary conditions to avoid the singular- 
ity. Physically speaking, this form (1) can be interpreted 
as follows. In the high velocity regime, where the radiative 
diffusion may become anisotropic in the comoving frame, 
the 'apparent' optical depth r would be of the order of 
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(2) 



That is, as the flow is accelerated and approaches to the 
speed of light, the optical depth becomes zero (outward 
peaking). In this case, the form (1) can be read as 

which recovers a similar form of a usual variable 
Eddington factor (see, e.g., Tamazawa et al. 1975). 
Hence, the applicability and accuracy of the form (1) from 
a low speed regime to a high speed one would be similar 
to those of a variable Eddington factor from an optically 
thick regime to an optically thin one. 

Similarly, for a spherically symmetric case, Akizuki and 
Fukue (2006) proposed a variable Eddington factor, which 
depends on the flow velocity (3 as well as the optical depth 
r: 



7(l + /3) + r 
7(l + /3) + 3r' 
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where 7 (= l/-\/l — P"^) is the Lorentz factor. 

In Fukue (2006b), the plane-parallel case was examined 
as an example, although the gravity of the central ob- 
ject was ignored for simplicity. In this paper, we thus 
consider the radiatively driven vertical outflow - moving 
photosphere - in a luminous flat disk within the frame- 
work of radiation transfer in the rclativistic regime using 
/(/3), while taking into account the gravity of the central 
object, although the gas pressure is ignored. 

In the next section we describe the basic equations in 
the vertical direction. In section 3 we show our numerical 
examination of the radiative flow. In section 4 we briefly 
apply the present model to the case of a supercritial ac- 
cretion disk. The final section is devoted to concluding 
remarks. 

2. Basic Equations and Boundary Conditions 

Let us suppose a luminous flat disk, inside of which 
gravitational or nuclear energy is released via viscous heat- 
ing or other processes. The radiation energy is trans- 
ported in the vertical direction, and the disk gas, itself, 
also moves in the vertical direction due to the action of 
radiation pressure. For simplicity, in the present paper, 
the radiation field is considered to be sufficiently intense 
that the gas pressure can be ignored: tenuous cold normal 
plasmas in the super-Eddington disk, cold pair plasmas 
in the sub-Eddington disk, or dusty plasmas in the sub- 
Eddington disk. As for the order of the flow velocity v, 
we consider the fully rclativistic regime, where the terms 
are retained up to the second order of (v/c). 



2.1. 



isic Equations 



Under these assumptions, the radiation hydrodynamic 
equations for steady vertical (z) flows without rotation are 
described as follows (Kato et al. 1998; Fukue 2006b). 

The continuity equation is 



pcu = J (= const.). 



(5) 



where p is the proper gas density, u the vertical four ve- 
locity, J the mass-loss rate per unit area, and c the speed 
of light. The four velocity u is related to the proper three 
velocity v hy u = jv/c, where 7 is the Lorentz factor, 

7 = yrT^=i/v/i-(w/c)2. 

The equation of motion is 
du GMz 



dz 



{R-r^fR 



[F7(l -f 2^2) - {cE + cP)-f'^u\ i&) 



where M is the mass of the central object, R ~ Vr^ + z^, 
r being the radius, (= 2GM/c^) is the Schwarzschild 
radius, Kabs and Kgca are the absorption and scattering 
opacities (gray), defined in the comoving frame, E the 
radiation energy density, F the radiative flux, and P the 
radiation pressure observed in the inertial frame. The 



first term in the square bracket on the right-hand side of 
equation (6) means the radiatively-driven force, which is 
modified to the order of v?' , whereas the second term is the 
radiation drag force, which is also modified, but roughly 
proportional to the velocity. As for the gravity, we adopt 
the pseudo-Newtonian potential (Paczyhski, Wiita 1980). 

When the gas pressure is ignored, the advection terms 
of the energy equation are dropped (cf. Kato et al. 1998), 
and heating is balanced with the radiative terms, 

= q+ - /9 (j - KabsC-E'7^ - KabsCPu^ + 2Kabs-F7") . (7) 

where is the internal heating and j is the emissivity 
defined in the comoving frame. In this equation (7), the 
third and fourth terms on the right-hand side appear in 
the relativistic regime. 

For radiation fields, the zeroth- moment equation be- 
comes 



dF r 



KahsCE + Ksca(c£^ + cP)u^ 

The first-moment equation is 



(8) 



dP 
dz 



PI 
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:a 

F{1 + 2m^) -I- KscnicE^ 



- cP)ju\ 



(9) 



In order to close moment equations for radiation fields, 
we need some closure relation. Instead of the usual 
Eddington approximation, we here adopt a velocity- 
dependent variable Eddington factor /(/3), 



Po = fW)Eo 



(10) 



in the comoving frame, where Pq and Eq are the quantities 
in the comoving frame. If we adopt this form (10) as the 
closure relation in the comoving frame, the transformed 
closure relation in the inertial frame is 

cP {l + u^- fu^) = cE {fj^ - u^) + 2F-fu (1 - /) , (11) 
or equivalently, 

cP (1 - fl3^) = cE (/ - /32) + 2FP (1 - /) . (12) 
As a form of the function /(/3) we adopt the simplest one: 



(13) 



for a plane-parallel geometry (Fukue 2006b; cf. Akizuki 
and Fukue 2006 for a spherically symmetric geometry). 

Eliminating j and E with the help of equations (7) and 
(11), equations (6), (8), and (9) become 



cp'u 



du 



c r„z 



2{R-r^YR 



Kabs + Ksca F{f-f^+U^)-cP{l + f)-1U 

7 



dF 
dz 



FUl^ + u')-cP{l + f)-fu 
/7^ — 



(14) 



, (15) 
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dP 
dz 



P : 7 — — — ■ (16) 



c — 
Introducing the optical depth by 
dr = -(^abs + Ksci>,)pdz, 



(17) 



and using continuity equation (5), equations (14)-(17) are 
rearranged as 



T du 
c J— = 



dr Kabs + Ksca2(i?-rg)2i? 

F{f^^ + u^)-cP{l + f)ju 



-7- 



J 



dz 



cu 



dr 
dF_ 
dr 



dP 



(/^abs H" ^sca)p 

F{fj^ + u^)-cP{l + f)^u 

^ -f 2 2 ' 

77 -u 
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('^abs + ^sca)p 

F{fj^ + u^)-cPil + f)ju 
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(18) 
(19) 

(20) 
(21) 



In this paper we assume that the heating takes place 
deep inside the disk at the midplane and in the atmosphere 
(7+ = 0. However, it is straightfoward to consider more 
general cases, where, e.g., the heating g"*" is proportional 
to the gas density p (cf. Fukue 2005a, b). 

We solved equations (18)-(21) for appropriate bound- 
ary conditions at the moving surface with a variable 
Eddington factor (13). 

2.2. Boundary Conditions 

As already pointed out in Fukue (2005a), the usual 
boundary conditions for the static atmosphere cannot be 
used for the present radiative flow, which moves with ve- 
locity at the order of the speed of light. 

When there is no motion in the gas ("static photo- 
sphere"), the radiation field just above the surface un- 
der the plane-parallel approximation is easily obtained. 
Namely, just above the disk with surface intensity Jg, the 
radiation energy density Eg , the radiative flux Fs , and the 
radiation pressure Pg are (2/c)7r/s, nig, and (2/3c)7r/s, 
respectively, where the subscript s denotes the values at 
the disk surface. However, the radiation field just above 
the surface changes when the gas itself does move upward 
("moving photosphere"), since the direciton and inten- 
sity of radiation change due to relativistic aberration and 
Doppler effect (cf. Kato et al. 1998; Fukue 2000). 

Let us suppose a situation that a flat infinite photo- 
sphere with surface intensity Is in the comoving frame is 
not static, but moving upward at a speed Vg (= c/?s, and 
the corresponding Lorentz factor is 7s). Then, just above 
the surface, the radiation energy density Eg, the radiative 
flux Fg, and the radiation pressure Pg measured in the 



inertial frame become, respectively. 



cK = 2ttL 



K - 2t:L 



cP, = 2ttL 
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(22) 
(23) 
(24) 



where Ug (= 7sVs/c) is the flow four velocity at the surface 
(Fukue 2005a). 

Thus, we impose the following boundary conditions: At 
the flow base (deep "inside" the atmosphere) with an ar- 
bitrary optical depth tq, the flow velocity u is zero, the 
radiative flux is Fq (which is a measure of the strength of 
radiation field), and the radiation pressure is Pq (which 
connects with the radiation pressure gradient and relates 
to the internal structure), where the subscript denotes 
the values at the flow base. At the flow top ("surface" 
of the atmosphere) where the optical depth is zero, the 
radiation field should satisfy the values above a moving 
photosphere given by equations (22)-(24): i.e.. 



Fg 



2 + 6pg + 6P^ 

3 + 8/3s + 3/32 



(25) 



where (3s is a final speed at the disk surface. 

Physically speaking, in the radiative flow starting from 
the flow base with an arbitrary optical depth tq , for initial 
values of Fq and Pq at the flow base, the final value of 
the flow velocity Vs at the flow top can be obtained by 
solving basic equations. Furthermore, the mass-loss rate J 
is determined as an eigenvalue so as to satisfy the bondary 
condition (23) at the flow top (cf. Fukue 2005a in the 
subrelativistic regime) . 

However, the permitted region for J is very tight, and it 
is difficult to search the value of J. Hence, in this paper, 
as a mathematically equivalent way, we fix the value of J, 
and search the value of Pa so as to satisfy the boundary 
condition (23). 

3. Relativistic Radiative Flow under Gravity 

In this section we show the relativistic radiative vertical 
flow in the luminous disk under the influence of gravity 
of the central object. In order to obtain the solution, 
as already stated, we numerically solve equations (18)- 
(21), starting from r = tq at 2: = with appropriate initial 
conditions for v, F, and P, down to r = so as to satisfy 
the boundary conditions (25) there. The parameters are 
the initial radius r on the disk, the initial optical depth 
To, which relates to the disk surface density, the initial 
flux Fq, which is the measure of the strength of radiation 
field to gravity, and the initial radiation pressure Pq at 
the disk base, which connects with the radiation pressure 
gradient in the vertical direction and relates to the disk 
internal structure. The mass-loss rate J is determined as 
an eigenvalue of the boundary condition at the fiow top. 

Several examples of numerical calculations arc shown in 
figures 1-3. Physical quantities arc normalized in terms of 
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Fig. 1. Flow velocity v (thick solid curve), flow height z 
(chain-dotted one), radiative flux F (dashed one), and radia- 
tion pressure P (dotted one), as a function of the optical depth 
T. The parameters are r = 3, tq = 1, Fq = 1, and Pq = 1.23. 
The mass- loss rate is J = 1. The quantities are normalized in 
units of c, rg, and LE/(47rr|). 

the speed of light c, the Schwarzschild radius rg, and the 

Eddington luminosity Le [= 47rcGM/(Kabs + '*sca)]. The 
units of F, cP, and c^J is LE/(47rrg). It should be noted 
that the solutions can be obtained for arbitrary optical 
depths To at the flow base (see Fukue 2006b, and section 
4). We here show, however, the cases of tq = 1, where the 
velocity change is remarkable. 

In figure 1 we show the flow velocity v (solid curve), 
the flow height z (chain-dotted one), the radiative flux F 
(dashed one), and the radiation pressure P (dotted one), 
as a function of the optical depth r for r = 3, tq^I, F^^l, 
and Po = 1-23 (i.e., J= 1). 

As the optical depth decreases from the flow base at 
the disk equator to the flow top at the disk surface, the 
radiative flux slightly decreases while the flow velocity in- 
creases; the radiative energy is converted to the flow bulk 
motion in the vertical direction. As usually known, in 
a static plane-parallel atomosphere, under the radiative 
equilibrium with gray approximation, the vertical flux F 
is conserved without heating source. In the present rela- 
tivistically moving atmosphere^ on the contrary, the radia- 
tive flux F decreases via Fu term, which acts to drive the 
gas toward the vertical direction. In the case of figure 1, 
the initial flux (Fq = 1) is nearly the local Eddington one, 
and therefore, the final fiow velocity is mildly relativistic 
due to the effect of gravity of the central object. Other 
parameter dependences are shown in figures 2 and 3. 

In figure 2 the flow velocity v are shown for several 
parameter set: A thick solid curve is for the typical case 
of r = 3, To = 1, Fo = 1, and Pq = 1.23 (i.e., J = 1). A 
solid curve is for the case of r = 3, tq = 1, Fq = 10, and 
Pq = 10.6 (i.e., J = 1), while a dashed one is for the case 
of r = 3, To = 1, Po = 1, and Pq = 1.041 (i.e., J = 0.005). 

As is easily expected, the flow velocity remarkably in- 
creases when the initial flux Pq is large (a solid curve). 
When the initial flux is small, on the other hand, the flow 
velocity becomes small. 

Even for the same initial flux, when the mass-loss rate 



Fig. 2. Flow velocity v as a function of the optical depth 
T. A thick solid curve is for the typical case of r = 3, tq = 1, 
Fo = 1, and Pq = 1.23 (i.e., J = 1). A solid curve is for the 
case of r = 3, tq = 1, Fq = 10, and Pq = 10.6 (i.e., J = 1), 
while a dashed one is for the case of r = 3, tq = 1, Po = 1> and 
Po = 1.041 (i.e., J = 0.005). The quantities are normalized in 
units of c, rg , and / (47rrg ) . 



2.0, 



1.0- 




Fig. 3. Flow velocity as a function of the optical depth t 
for several initial radii r: r = 3 (a thick solid curve), r = 2 (a 
solid one), and r = 1.5 (a dashed one). Other parameters are 
TO = 1, Po = 1, and Po = 1.23 (i.e., J = 1). The quantities are 
normalized in units of c, rg, and LE/(47rrg). 

is small, the flow velocity remarkably increases (a dashed 
curve). This is because for a small mass- loss rate the 
density decreases, and therefore, the vertical height z be- 
comes large, compared with the case for a large mass-loss 
rate with the same optical depth. As a result, the gas 
is accelerated along the long distance to reach the highly 
relativistic regime. 

In flgure 3 the dependence on the initial radius r is 
shown: r = 3 (a thick solid curve), r = 2 (a solid one), 
and r = 1.5 (a dashed one). Other parameters are tq = 1, 
Po = 1, and Po = 1.23 (i.e., J = 1). 

For the same initial conditions, the flow velocity de- 
creases as the initial radius decreases. This is just the ef- 
fect of gravity of the central object. In the case of r = 1.5, 
the velocity is slightly decelarated toward the surface due 
to gravity. 
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4. Relativistic Radiative Flow in the Critical 
Disk 

In this section we apply the present model to the mass 
outflow in the luminous supercritical accretion disks (cf. 
Fukue 2006a for a subrelativistic case) . 

When the mass-accretion rate M in the disk around 
a central object of mass M highly exceeds the critical 
rate Merit, defined by Merit = -^e/c^, the disk is believed 
to be in the supercritial regime, and the disk luminos- 
ity exceeds the Eddington one. Such a supercritical ac- 
cretion disk, a so-called slim disk, has been extensively 
studied, both numerically and analytically (Abramowicz 
et al. 1988; Eggum et al. 1988; Szuszkiewicz et al. 1996; 
Bcloborodov 1998; Watarai, Fukue 1999; Watarai et al. 
2000; Mineshige et al. 2000; Fukue 2000; Kitabatake et 
al. 2002; Okuda 2002; Ohsuga et al. 2002, 2003, 2005; 
Watarai, Mineshige 2003; Fukue 2004; Ohsuga 2006). It 
was found that the optically-thick supercritical disk is 
roughly expressed by a self-similar model (e.g., Watarai, 
Fukue 1999; Fukue 2000; Kitabatake et al. 2002; Fukue 
2004). Except for the case of Fukue (2004), however, many 
of these analytical models did not consider the mass out- 
flow from the disk surface. Hence, in this paper we adopt 
the model developted by Fukue (2004), as a background 
supercritical disk model. 

In the critical accretion disk constructed by Fukue 
(2004), the mass-accretion rate is assumed to be regulated 
just at the critical rate with the help of wind mass-loss. 
Outside some critical radius, the disk is in a radiation- 
pressure dominated standard state, while inside the criti- 
cal radius the disk is in a critical state, where the excess 
mass is expelled by wind and the accretion rate is kept to 
be just at the critical rate at any radius. Here, the critical 
radius is derived as 

9'\/3crT 



-Mi 



input ' 



(26) 



where Mi„put is the accretion rate at the outer edge of the 
disk, and m — Minput/Mcrit- Outside rer, the accretion 
rate is constant, while, inside rer the accretion rate would 
vary as 



M(r) 



IQlTCm-n 



■ Mi 



input ~ 



(27) 



In such a critical accretion disk, the disk thickness H is 
conical as 

H ^ If m 

— = = - In IH 

r ^ ^ 4 V 20 



(28) 



where C3 is some numerical coefficient determined by the 
similar procedure in Narayan and Yi (1994) for optically- 
thin advection-dominated disks. The second equality 
of this relation comes from the numerical calculation 
(Watarai et al. 2000). Although the mass loss was not 
considered in Watarai et al. (2000), we adopted this re- 
lation as some measure: when the normalized accretion 
rate m is 1000, the coefficient y/c^ becomes 0.983. 
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Fig. 4. Flow velocity v as a function of r for several values 
of r. The values of r are, from left to right, 3, 4, 5, 6, and 
7. The quantities are normalized in units of c and rg. The 
parameters of the critical disk is m = 1000 and « = 1. 



Furthermore, in Fukue (2004), several alternatives are 
discussed, and some of them gives the physical quantities 
of the critical accretion disk with mass loss as 



16V6 




(29) 
(30) 
(31) 



where a is the viscous parameter. 

In the present non-dimensional unit in terms of c, rg, 
and the Eddington luminosity Le, these physical quanti- 
ties are expressed as 



H^^^r 
16^/6 



To 

Po = 



a 
1V6 

a 



,3/2 ■ 



(32) 
(33) 

(34) 

(35) 



where the symbol "hat" (say, f) is dropped. 

Using these relations, we can solve the basic equations, 
and obtain numerical solutions at each radius r. The ex- 
ample in the case of rh = 1000 and a = 1 is shown in figures 
4 and 5. 

In figure 4 we show the flow velocity w as a function of 
optical depth r for several values of r. The quantities are 
normalized in units of c and rg. The parameters of the 
critical disk is m = 1000 and a — 1. 

As can be seen in figure 4, the flow velocity v varies self- 
similarly for different values of radii. This may be reflected 
the initial self-similar models. As a result, in each radius 
with different optical depth, the flow final speed is almost 
same. 

In figure 5 we show several quantities for each radius r: 
The disk height H (dashed curve) and the optical depth 
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Fig. 5. Several quantities for eacii radius r: The disk heiglit 
H (dashed curve) and the optical depth r (chain-dotted one) 
are from the critical model, while the height Zs and velocity Vs 
(solid curves) at the flow top and the mass-loss rate J (dashed 
one) are the results of the present numerical calculations. The 
quantities are normalized in units of c and rg. The parameters 
of the critical disk is rri = 1000 and a = l. 

To (chain-dotted one) are from the critical accretion disk 
model, while the height Zg and velocity Vg (solid curves) 
at the flow top and the mass-loss rate J (dashed one) are 
the results of numerical calculations. The quantities are 
normalized in units of c and rg. The parameters of the 
critical disk is m = 1000 and a = l. 

As can be seen in figure 5, the flow height Zg is enor- 
mously large. Hence, rigorously speaking, the present 
plane-parallel approximation violates in this application 
of m = 1000, and two-dimensional numerical simulation 
should be needed in such a case. However, we can see 
several insights from the present case. 

First, the final speed of the fiow accelerated in the lu- 
minous critical disk becomes sufficiently relativistic. In 
other words, relativistic jets can form in such a luminous 
accretion disk. In addition, this final speed does not de- 
pend on the initial radius so much due maybe to the initial 
self-similarity. Second, on contrary to the final speed, the 
mass-loss rate per unit area increases as the initial radius 
decreases; it is roughly approximated by 6/r. On the 
other hand, the model mass-loss rate (Fukue 2004. 2006a) 
becomes J — 1/r, that is qualitatively same, but quantita- 
tively different from the present numerical values. Hence, 
the mass loss from the critical disk may be concentrated in 
the inner region, although the true mass-loss rate cannot 
be determined at the present simple state. This nature, 
however, is also convenient for centrally concentrated jets. 

5. Concluding Remarks 

In this paper we have examined the relativistic radia- 
tive transfer fiow in a luminous disk in the relativistic 
regime of (w/c)^, taking account of gravity of the cen- 
tral object. In such a relativistic regime, we adopt the 
velocity-dependent variable Eddington factor. The fiow is 
assumed to be vertical, and the gas pressure is ignored for 
simplicity. The basic equations are numerically solved as 
a function of the optical depth r, and the flow velocity v, 



the height z, the radiative flux F, and the radiation pres- 
sure P are obtained for a given radius the initial optical 
depth To, and the initial coditions at the fiow base (disk 
"inside" ), whereas the mass- loss rate J is determined as an 
eigenvalue of the boundary condition at the fiow top (disk 
"surface"). For sufficiently luminous cases, the flow re- 
sembles the case without gravity. For less-luminous cases, 
however, the flow velocity decreases. 

Application to the critical accretion disk was also exam- 
ined. If the disk thickness becomes so large, the present 
plane-parallel approximation violates and other treat- 
ment, such as numerical simulations, should be needed. 

The radiative flow investigaed in the present paper must 
be a quite fundamental problem for accretion-disk physics 
and astrophysical jet formation, although there are many 
simplifications at the present stage. For example, we have 
ignored the gas pressure. In general cases, where the gas 
pressure is considered, there usually appears sonic points, 
and the flow is accelerated from subsonic to supersonic. 
In this paper we consider a purely vertical flow, and the 
cross section of the flow is constant. If the cross section of 
the flow increases along the flow, the flow properties such 
as a transonic nature would be influenced. Moreover, we 
do not consider the rotation of the gas. In accretion disks 
around a black hole, the gas usually rotates around the 
hole at the relativistic speed. In the vicinity of the equa- 
tor, the vertical flow approximation safely holds, since the 
radial gravity is balanced with the centrifugal force. When 
the fiow is accelerated to be lift up the large height, the 
streamline would be curved outward, since the centrifu- 
gal force dominates. This would violate the vetical fiow 
approximation. 

There remain many problems to be solved. 
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